BOUNDARY CHARACTERISTIC POINT REGULARITY FOR 
SEMILINEAR REACTION-DIFFUSION EQUATIONS: 
TOWARDS AN ODE CRITERION 



V.A. GALAKTIONOV AND V. MAZ'YA 

Abstract. The classic problem of regularity of boundary characteristic points for semi- 
linear heat equations with homogeneous Dirichlet conditions is considered. It is shown 
that famous Petrovskii's (the so-called, 2-y/log log) criterion of boundary regularity for 
the heat equation (1934) can be adapted to classes of semilinear parabolic equations 
of reaction-diffusion type, and now this takes the form of an ODE regularity criterion. 
Namely, after a special matching with a Boundary Layer, the regularity problem reduces 
to a one-dimensional perturbed nonlinear dynamical system for the first Fourier-like 
coefficient of the solution in an Inner Region. 

A similar ODE criterion, with an analogous matching procedures, is shown formally 
to exist for semilinear fourth-order bi-harmonic equations of reaction-diffusion type. Ex- 
tensions to regularity problems of backward paraboloid vertexes in M. N are discussed. 



1. Introduction: main semilinear equations, classic regularity problem, 

results, and layout 

1.1. Semilinear reaction-diffusion PDEs near parabola vertexes. The present pa- 
per is devoted to a systematic study of the regularity of the origin (0,0) as a boundary 
point for semilinear heat (reaction-diffusion) equations, which we first consider in ID: 

(1.1) u t = u xx + f(x,t,u) in Q C R x [-1,0), u = on dQ , u(x, -1) = u Q (x). 



Here Q is a sufficiently smooth domain such that (0, 0) G dQ (dQ denotes the lateral 
boundary of Q ) is its only characteristic boundary point, i.e., in the {x, £}-plane, 

(i) the straight line {t = 0} is tangent to dQ at this point, and 

(ii) no such points exist on dQo for t G [—1,0). 

According to (11.11) . we pose the zero Dirichlet condition on the lateral boundary 8Qq and 
prescribe arbitrary bounded initial data uq(x) at t = — 1 in Q fl {t = — 1}. 

We naturally assume that nonlinearities f(x,t,u) in (II. ip satisfy necessary regularity 
and growth in u hypotheses that guarantee existence and uniqueness of a smooth classical 
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solution u G Cf'^(Qo) H C(Qo)0 of ( 11. ip in Qo by classic parabolic theory; see e.g., well- 
known monographs [U [121 HI HZ]- A standard regularity issue in general PDE theory is 
then as follows: how to control the solution at the vertex, i.e., the problem is: 

u(0,(T) =? 

Note that this does not a priori exclude blow-up at the vertex (regardless zero Dirichlet 
conditions on the lateral boundary) where |m(0,0~)| = +00, in the sense of limsup. 

More precisely, as customary in regularity theory, the goal is to derive conditions show- 
ing how given smooth nonlinear perturbations /(•) can affect regularity of the vertex (0, 0) 
of a such a backward parabola dQo- The regularity of (0,0) (in Wiener's sense) means: 

(1.2) regularity of (0,0) : -u(0,0 _ ) = for any initial data ito- 

Recall that, as is well known, for nonlinearities /(•) = 0, i.e., for the pure heat equation 

(1.3) u t = u xx in Q , u = on dQ , u(x, -1) = u (x), 

this regularity problem was solved by Petrovskii in 1934 [53] 154"] . who introduced his 
famous Petrovskii 's regularity criterion (the so-called "2^/log log-one"; see details below). 

Indeed, many and often strong and delicate boundary regularity and related asymp- 
totic results are now known for a number of quasilinear parabolic equations, including 
even a few for degenerate porous medium operators. Nevertheless, some difficult ques- 
tions remain open even for the second-order parabolic equations with order-preserving 
semigroups. We refer to the results and surveys in [Tj [2j [3j HH [27J EJ ETJ as a guide 
to a full history and already existing interesting extensions of these important results. 
Concerning further developing of Wiener's ideas in linear parabolic equations, see refer- 
ences and results in [UJJ |39j [TJ] and in [57]. However, a more systematic study of those 
regularity issues for equations such as (II. ip with rather general nonlinear perturbations 
/(•) was not done properly still. In fact, it turned out that, for such arbitrary /(-)' s > the 
classical barrier methods hardly applied and another asymptotic approach was necessary. 
We propose this in the present paper for a wide class of semilinear parabolic PDEs. 

1.2. Layout of the paper: key models, nonlinearities, and extensions. Section 
[2] contain some preliminary discussions and results. In Sections [3HS1 the main goal is to 
show how a general "nonlinear perturbation" /(■) in (II. ip affects the regularity condi- 
tions by deriving sharp asymptotics of solutions near characteristic points. To this end, 
we apply a method of a matched asymptotic (blow-up) expansion, where the Boundary 
Layer behaviour close to the lateral boundary dQ (Section [1J is matched, as t — > 0~, 
with a centre subspace behaviour in an Inner Region (Section [5]). This leads to a non- 
linear dynamical system for the first Fourier coefficient in the eigenfunction expansion 
via standard Hermite polynomials as eigenf unctions of the linear Hermite operator ob- 
tained after blow-up scaling near the vertex. Overall, the vertex regularity is shown to 

1 As customary in parabolic PDEs [16] , the closure Qo does not include the "upper lid" , which is the 
vertex (0, 0) only. 
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be governed by an ODE criterion, which principally does not admit any simply integral 
(Osgood-Dini-type) treatment as in Petrovskii's one. 

Indeed, such an approach falls into the scope of typical ideas of asymptotic PDE theory, 
which got a full mathematical justification for many problems of interest. In particular, 
we refer to a recent general asymptotic analysis performed in [31]. According to its clas- 
sification, our matched blow-up approach corresponds to perturbed ONE-DIMENSIONAL 
dynamical systems, i.e., to a rather elementary case being however a constructive one that 
detects a number of new asymptotic/regularity results. 

In particular, to show a typical "interaction" between the linear Laplacian and the 
nonlinear perturbations in (11. ip . we initially concentrate on the simplest case, with 

(1.4) /(•)= (Zj)«(«)« for te[-l,0), 
where k(u) is a smooth enough function satisfying 

(1.5) k(u) — > as u — > 0, \k(u)\ < 1, k(u) ^ for u^O. 

We show that the nonlinear perturbation (II. 4p will then affect Petrovskii's 2i/log log- 
time-factor starting from some awkward looking functions such as 

(1.6) k(u) ~ |ln|u| |!e-( 3 v^l ln HI) 2/3 -> as u -> 0. 

For more general nonlinearities, we derive a so-called ODE regularity criterion of the 
vertex (0,0), meaning that a special nonlinear ODE for the first Fourier coefficients of 
rescaled solutions takes responsibility for the vertex regularity/irregularity. 

The present research has been inspired by the regularity study of quasilinear ellip- 
tic equations with quadratic gradient-dependent nonlinearities [4"S] . where, in 2D, new 
asymptotics of solutions near corner points were discovered. We also refer to monographs 
[Ml ESI EH SSI HZ] and [33], [12] [II] as an update guide to elliptic regularity theory 
including higher-order equations. Sharp asymptotics of solutions of the heat equation 
in domains with conical points were derived in [251 E2l ESI [31]. Higher-order parabolic 
equations were treated in [30, 29J. It turned out that, unlike the present study, such 
asymptotics are of a self-similar form. See also [38J for a good short survey including 
compressible/incompressible Stokes and Navier-Stokes problems. 

Therefore, as a next key model regularity problem, we briefly reflect the main differences 
and difficulties, which occur by studying the regularity issues for parabolic equations with 
a typical quadratic gradient dependence in the nonlinear term: 

(1.7) u t = u xx + k{u)u {u x ) 2 in Q , u = on dQ , u(x, — 1) = u Q (x). 

Then the ODE regularity criterion is expressed in terms of another ID dynamical system, 
with a weaker nonlinearity. We then convincingly show that, for any k(u) in (II. 7p . sat- 
isfying (II. 5p . Petrovskii's linear regularity criterion takes place, i.e., remains the same as 
for the heat equation (II. 3p . 

We also pay some attention to extensions to similar regularity problems in domains 
Qo C x [—1, 1), with dQo having a backward paraboloid shape and the vertex (0,0) 
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being their characteristic point. In Section |3l we thus discuss the semilinear problems: 

in Q , u = on dQ . 

n(u)u\Vu\ 2 

Finally, in Appendix B (Appendix A is devoted to the corresponding spectral theory of 
rescaled operators), we show how our approach can be extended to higher-order PDEs, 
e.g., for the semilinear bi-harmonic equations having similar nonlinearities, with also 
zero Dirichlet conditions on dQo and bounded initial data u in Q fl {t = —1}. The 
mathematical analysis becomes much more difficult and we do not justify rigorously all its 
main steps such as the boundary layer and matching with the Inner Region asymptotics. 
Moreover, the ID dynamical system for the first Fourier coefficients becomes also more 
delicate and does not admit such a complete analysis, though some definite conclusions 
are possible. We must admit that this part of our study is formal, though some steps are 
expected to admit a full justification, which nevertheless can be rather time-consuming. 



2. Petrovskii's 2 A/log log— criterion of 1934 and some extensions 

We need to explain some details of Petrovskii's classic regularity analysis for the heat 
equation performed in 1934-35. Following his study, we consider the one-dimensional 
case N — 1, where the analysis becomes more clear. Moreover, our further extensions 
to bi-harmonic operators (Appendix B) will be also performed for N = 1, in view of 
rather complicated asymptotics occurred, so we are not interested in involving extra 
technicalities. 

After Wiener's pioneering regularity criterion for the Laplace equation in 1924 [58J, 
I.G. Petrovskii [531 Ell was the nrs ^ wno completed the study of the regularity question 
for the ID and 2D heat equation in a non-cylindrical domain. We formulate his result in 
a blow-up manner, which in fact was already used by Petrovskii in 1934 [53J. 

He considered the question on an irregular or regular vertex (x, t) = (0, 0) in the initial- 
boundary value problem (IBVP) 

{u t = u xx in Q = {\x\ < R(t), -1 < t < 0}, R(t) -> + as t ->■ 0", 
with bounded smooth data u(x,0) = Uq(x) on [—R{—i),R(—l)\. 

Here the lateral boundary {x = ±R(t), t e [—1,0)} is given by a function R(t) that is 
assumed to be positive, strictly monotone, C 1 -smooth for all —1 < t < (with R'(t) > 0), 
and is allowed to have a singularity of R'(t) at t = 0~ only. The regularity analysis then 
detects the value of u(x, t) at the end "blow-up" characteristic point (0, 0~), to which the 
domain Q "shrinks" as t — > 0~ . 

Remark: on first parabolic regularity results for m — 1 and m > 2. It is well- 
known that, for the heat equation, the first existence of a classical solution (i.e., continuous 
at (0, 0)) was obtained by Gevrey in 1913-14 [23] (see Petrovskii's references in [531 P- 55] 
and [541 p. 425]), which assumed that the Holder exponent of R(t) is larger than |. In 
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our setting, at t — , this comprises all types of boundaries given by the functions: 

(2.2) R(t) = {-t) v with any v > \ are regular (Gevrey, 1913-14). 

For 2mth-order parabolic poly-harmonic equations such as (IA.3j) below, a similar result: 

(2.3) for R(t) = (— t) 5 ™, the problem is uniquely solvable 

was proved^ by Mihailov [1SJ almost sixty years later and fifty years ago. 
Definitions, (i) Regular point: as usual in potential theory, the point (x, t) = (0, 0) 
is called regular (in Wiener's sense, see [13]), if any value of the solution u(x,t) can be 
prescribed there by continuity as a standard boundary value on 8Qq. In particular, as a 
convenient and key for us evolution illustration, (0,0) is regular if the continuity holds for 
any initial data u (x) in the following sense: 

(2.4) u = at the lateral boundary {\x\ = R(t), -1 < t < 0} u(0, (T) = 0. 

(ii) Irregular point: otherwise, the point (0,0) is irregular, if the value u(0,0~) is not 
fixed by boundary conditions, i.e., w(0,0) ^ for some data u , and hence is given by a 
"blow-up evolution" as t — > 0~. Hence, formally, (0,0) does not belong to the parabolic 
boundary of Q . 



Petrovskii's '^^log log". Using novel barriers as upper and lower solutions of (12. ip . 
Petrovskii [521 El] established the following "2-\/log log-criterion" : 

j)R(t) = 2 v /3 t v/ln|ln(-£)| (0,0) is regular, and 

'ii) R(t) = 2(1 + e)^i ^ln | ln(-t)|, e > =^ (0, 0) is irregular. 



More precisely, he also showed that, for the curve expressed in terms of a positive 
function p(h) — > + as h — > + (p(h) ~ t^-^j is about right) as follows: 



(2.6) R(t) = v/-lnp(-t), 

the sharp regularity criterion holds (in Petrovskii's original notation): 



(2.7) 



j P (h)^\in P (h)\ dh< ^ + OQ ^ (Q) -) ig irregular ( re g U i ar ). 



Both converging (irregularity) and diverging (regularity) integrals in (12. 7\\ as Dini- 
Osgood-type regularity criteria already appeared in the first Petrovskii paper [531 P- 56] 
of 1934. Further historical and mathematical comments concerning Petrovskii's analysis 
including earlier [28] (1933) Khinchin's criterion in a probability representation can be 
found in a survey in [IS] . 



2 However, in Slobodetskii-Sobolev classes, i.e., continuity at (0,0) is not understood in the above 
Wiener classic sense. In fact, for m = 2, Wiener's one (|2.4p fails for the parabola with R(t) = 5(—t)* in 
, while the one with R(t) = 4(— t)z remains regular, [TSJ § 4.3]. 
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Petrovskii's integral criterion of the Dini-Osgood type given in ( 12 .7p is true in the 
N- dimensional radial case with (see [TJ [2J [3] for a more recent updating) 

I K 

(2.8) a/ | \np(h)\ replaced by | lnp(h)\ 2 . 

It is worth mentioning that, as far as we know, (12. 5p is the first clear appearance of the 
"magic" y/log log in PDE theory, currently associated with the "blow-up behaviour" of 
the domain Q and corresponding solutions. Concerning other classes of nonlinear PDEs 
generating blow-up y/log log in other settings, see references in [T5]. 

Thus, since the 1930s, Petrovskii's regularity ^log log-factor entered parabolic theory 
and generated new types of asymptotic blow-up problems, which have been solved for a 
wide class of parabolic equations with variable coefficients as well as for some quasilinear 
ones. Nevertheless, such asymptotic problems were very delicate and some of them of 
Petrovskii's type remained open even in the second-order case, i.e., for (II. ip . to be solved 
in the present paper for the first time. 

3. Preliminaries of matched asymptotic expansion 

3.1. The basic initial-boundary value problem. Thus, we consider the semilinear 
parabolic equation (II. ip . with a simple, "basic" nonlinear perturbation, which we take in 
the separable form (II. 4p . The eventual ODE regularity criterion will then also include 
the behaviour of the nonlinear coefficient k(u) as u — > 0. 

Hence, our basic second- order initial-boundary value problem (IBVP) takes the form: 

{u t = u xx + k{u)u in Q = {\x\ < R(t), -1 < t < 0}, 
u = at x = ±R(t), — 1 < t < 0, 
u(x, 0) = u (x) on [-R(-1),R(-1)], 

where uo(x) is a bounded and smooth function, w (±i?(— 1)) = 0. We then apply to (13.11) 
both Definitions from the previous Sectional 

3.2. Slow growing factor </>(r). According to (12.51) . we need to assume that 

(3.2) R(t) = (-t)*(p(r), where r = - ln(-t) -»• +oo as t -»■ 0~. 
Here, v 9 ( r ) > is a smooth monotone increasing function satisfying <f'(r) > 0, 

(3.3) <p(r) -> +oo, <pf(r) -)> + , and ^ -> as r -> +oo. 

Moreover, as a sharper characterization of the above class of slow growing functions, we 
use the following criterion: 

(3.4) (0^)'^oo as r^+oo fo/(r) ^ 0). 

This is a typical condition in blow-up analysis, which distinguishes classes of exponential 
(the limit in ( 13. 4 p is 0), power-like (a constant ^ 0), and slow-growing functions. See 
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[55| pp. 390-400], where in Lemma 1 on p. 400, extra properties of slow-growing functions 
f)3.4p are proved. For instance, one can use a comparison of such a <£>(t) with any power: 

(3.5) for any a > 0, ^(t) < r" and (p'(r) < r a ~ l for r > 1. 

Such estimates are useful in evaluating perturbation terms in the rescaled equations. 

Thus, the monotone positive function <p(r) in ( 13. 2 p is assumed to determine a sharp 
behaviour of the boundary of Qo near the shrinking point (0, 0) to guarantee its regularity. 
In Petrovskii's criterion ( 12. 5p . the almost optimal function, satisfying ( 13. 3p . ( 13. 4ft . is 

(3.6) f*( T ) = 2Vln r as r — > +oo. 

3.3. First kernel scaling and two region expansion. By ( 13. 2p . we perform the sim- 
ilarity scaling 

(3.7) u(x,t) = v(y,r), where y = jz^72- 
Then the rescaled function v(y,r) now solves the rescaled IB VP 

{v T — B*f + k(v)v = v yy — | yv y + k(v)v in Q = {\y\ < <p(r), t > 0}, 
v = at y = ±^(r), r > 0, 
v(0,y) = v (y)=u (y) on [-i2(-l), i2(-l)]. 
The rescaled equation in (13.81) . for the first time, shows how the classic Hermite operator 
(3.9) B* = Dl - | yD y 

occurs after blow-up scaling (13. 7p . By the divergence (13.31) of <p(r) — > +oo as r — > 
+oo, it follows that sharp asymptotics of solutions will essentially depend on the spectral 
properties of the linear operator B* on the whole line M (see Appendix A), as well as 
on the nonlinearity k(v)v, so that such an asymptotic "interaction" between linear and 
nonlinear operators therein eventually determines regularity of the vertex. 

Studying asymptotics for the rescaled problem (13.81) . as usual in asymptotic analysis, 
this blow-up problem is solved by matching of expansions in two regions: 

(i) In an Inner Region, which includes arbitrary compact subsets in y containing the 
origin y — 0, and 

(ii) In a Boundary Region close to the boundaries y = ±y?(r), where a boundary layer 
occurs. 

Actually, such a two-region structure, with the asymptotics specified below, defines the 
class of generic solutions under consideration. We begin with the simpler analysis in the 
Boundary Region (ii). 

4. Boundary layer (BL) theory 

4.1. BL-scaling and a perturbed parabolic equation. Sufficiently close to the lateral 
boundary of Qo, it is natural to introduce the variables 

(4.1) z = and v(y, r) = w(z, r) w T = ^ w zz - \ zw z + ^ zw z + k(w)w. 
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We next introduce the BL-variables 



(4.2) f = <p 2 (r)(l - z) = tp(ip -y), (^ 2 (r)dr = ds, and w(z, r) = p(s)g(£, s), 

where p(s) > for s ^> 1 is an unknown scaling time-factor depending on the function 
(p(r). As usual, this p-scaling is chosen to get uniformly bounded rescaled solutions, i.e., 
for nonnegative solutions, 

(4.3) sup ? g(£, s) = 1 for all s > 1 

(for solutions which remain of changing sign for s ^> 1, one takes s)| in (14.3j) ) . By 
the Strong Maximum Principle (Sturm's Theorem on zero sets, see [5]), v y (y,r) has a 
finite number of zeros in y for any r > (possible supremum points) and a standard 
argument ensures that the normalization (14. 3 j) implies that such a p(s) can be treated as 
sufficiently smooth for s ^> 10. This describes the class of solutions under consideration. 
For instance, by the Maximum Principle, it is particular easier to work out, when: 

(4.4) (14.31) holds for all nonnegative solutions u(x, t) ^ 0. 

Respectively, for non-positive solutions, one can use —1 as the normalization in (14.31) . 

On substitution into the PDE in (14.11) . we obtain the following small nonlinear pertur- 
bation of a linear uniformly parabolic equation: 

(4 9s = A<? - | ^ £<?e - f (1 - 

-2^^-fg+^K(pg)g, where Ag = g" + \g>. 

As usual in boundary layer theory, this means that we then are looking for a generic 
pattern of the behaviour described by ( 14. 5 p on compact subsets near the lateral boundary, 

(4.6) |f | = o((p~ 2 {t)) => \z- 1| = o(ip- 4 (r)) as r -»■ +oo. 

On these space-time compact subsets, the second term on the right-hand side of ( 14. 5 p 
becomes asymptotically small, while all the other linear ones are much smaller in view of 
the slow growth/decay assumptions such as (13 ,4p for ip(r) and p(s). 

4.2. Passing to the limit and convergence to a BL-profile. Thus, we arrive at a 
uniformly parabolic equation ( 14. 5 P perturbed by a number of linear and nonlinear terms 
being, under given hypothesis, asymptotically small perturbations of the stationary elliptic 
operator A. In particular, the last nonlinear term in ( 14. 5 p is clearly asymptotically small 
by the hypotheses (I1.5P and ( 13. 3p . so that, for uniformly bounded g, 

(4.7) -4^y g n(p(s)g) -> as r -»■ +oo. 



3 On the other hand, one can normalize in (|4.2[) by the smooth function p(s) = v(0, r), which also can 
be regarded as positive (negative) for r>l (infinitely many sign changes of v(0, r) for t>1 also mean 
that v(y,r) has infinitely many zeros in y that is impossible for the heat equation [5]). This leads to 
some slight technical differences, though makes the normalization (|4.8p below more straightforward. 
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The BL representation (14. 2p . by using the rescaling and (14. 3p . naturally leads to the 
following asymptotic behaviour at infinity: 

(4.8) lim s _> +00 g(£, s) -)■ 1 as £ ->■ +oo, 

where all the derivatives also vanish. Then, we arrive at the problem of passing to the 
limit as s — > +oo in the problem (14. 5p . (14. 8p . Since, by the definition in (14. 2p . the rescaled 
orbit {g(s), s > 0} is uniformly bounded, by classic parabolic interior regularity theory 
[T7J [121 H3]> one can pass to the limit in (14. 5 p along a subsequence {sk} — > +oo. Namely, 
we have that, uniformly on compact subsets defined in (14. 6p . as k — > oo, 

(4.9) g(s k + s) — > h(s), where h s = Ah, h = at £ = 0, /i|^=+oo = 1- 
Consider this /zmzt (at s = +oo) equation obtained from (14. 5ft : 

(4.10) /i s = Ah = % + |^ in l + xl +i h(0, s) = 0, /i(+oo,s) = l. 

It is a linear parabolic PDE in the unbounded domain M + , governed by the operator A 
admitting a standard symmetric representation in a weighted space. Namely, we have: 

Proposition 4.1. (i) (I4.10p is a gradient system in a weighted L 2 -space, and 

(ii) for bounded orbits, the u -limit set Q of (I4.10p consists of a unique stationary profile 

(4.11) flb(0 = l-e-* /a , 

and flo is uniformly stable in the Lyapunov sense in a weighted L 2 -space. 
Proof. As a 2nd-order equation, (I4.10p is written in a symmetric form, 

(4.12) e^ 2 h s = (e^ 2 h^, 

and hence admits multiplication by h s in L 2 that yields a monotone Lyapunov function: 

(4.13) IhS^WK = -fe^ 2 (h s ) 2 dt < 0. 

Note that, in (I4.12p . the derivatives h^ and h s have to have an exponential decay at infinity 
in order the seminorms involved to make sense. It is essential that the limit profile (14. lip 
perfectly suits both. 

Thus, the problem (I4.5P is a perturbed gradient system, that allows to pass to the limit 
s — > +oo by using power tools of gradient system theory; see e.g., Hale [26J. 

(ii) For a given bounded orbit {h(s)}, denote h(s) = go + w(s), so that w(s) solves the 
same equation (I4.12p . Multiplying by w(s) in L 2 yields 

(4.14) \ A / 4l 2 w 2 dC = - / e«/ 2 (^) 2 de < 

for any nontrivial solutions, whence the uniform stability (contractivity) property. □ 

Finally, we state the main stabilization result in the boundary layer, which establishes 
the actual class of generic solutions we are dealing with. 



Theorem 4.1. (i) There exists a class of solutions of the perturbed equation (j4.5j) . for 

which, in a weighted L 2 -space and uniformly on compact subsets, 

(4.15) g(t,s)^-g (0 as s -> +00. 

(ii) ( I4.15P is particularly true for all nontrivial nonnegative solutions. 

Proof, (i) Under given hypotheses, the uniform stability result in (ii) of Proposition 14.11 
implies [22], Ch. 1] that the w-limit set of the asymptotically perturbed equation (14. 5 j) is 
contained in that for the limit one (14.101) . which consists of the unique profile (14. lip . 

(ii) This follows from the construction, since then p(s) in (14.21) can be chosen always 
positive. Then in the limit we are guaranteed to arrive at the gradient problem (14.91) 
admitting the unique uniformly stable stationary point (14.111) . □ 



5. Inner Region expansion: towards an ODE regularity criterion 

5.1. The Cauchy problem setting, eigenfunction expansion, and matching. In 

Inner Region, we deal with the original rescaled problem (I3.8p . Without loss of generality, 
again for simplicity of final, rather technical and involved calculations, we consider even 
solutions defined for y > by assuming the symmetry condition at the origin 

(5.1) v y = at y = 0. 

As customary in classic PDE and potential theory (see e.g., Vladimirov [23 § 6]), we 
extend v (y, r) by beyond the boundary points, i.e., for y > ip(r): 



(5.2) v(y, t) = v(y,r)H(<f(r) - y) 



v(y,r) for < y < y?(r), 
for y > <p(r), 



where H is the Heaviside function. Since v = on the lateral boundary {y = ip{r)}, one 
can check that, in the sense of distributions, 

(5.3) v T = v T H, v y = v y H, and v yy = v yy H - v y \ 5(y - ip). 
Therefore, v satisfies the Cauchy problem: 

(5.4) v T = B*v + v y \ y=v{T) 5(y - p(r)) + k(v)v in R x E+. 

Since, by construction, the extended solution (I5.2p is uniformly bounded in Zy,(R), 
we can use the converging in the mean (and uniformly on compact subsets in y) the 
eigenfunction expansion via the standard Hermite polynomials given in (1A.22I) for m = 1: 

(5-5) v(y, r) = E( fe >o) ak(r)i>* k {y). 

Actually, as follows from BL-theory from Section m Theorem 14.11 that the only possible 
solutions admitting matching with (14.151) possess a constant in y behaviour on compact 
subsets in y, i.e., 

(5.6) v(y, t) = (Zo(t) • 1(1 + o(l)) as r — > +00. 
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Indeed, this "1" well corresponds to the first Hermite polynomial ipo(y) = 1 in (15. 5p . 
Since Ao = for this "polynomial", the behaviour (I5.6P can be referred as to a "center 
subspace" one for the operator B* in (13.91) . though we do not use this fact at all. 

Thus, by BL-theory establishing the boundary behaviour (14.21) for r 3> 1, which we 
state again: in the rescaled sense, on the given compact subsets, 

(5-7) v(y, r) = p(s)g (^(r)(l - ^)) (1 + o(l)). 

Overall, in the class of generic solutions satisfying the BL-expansion, we concentrate on 
the first Fourier pattern associated with 

(5.8) k = 0: A = and = 1 (ip {y) = F(y), the Gaussian §KJM). 
The corresponding normalization condition is key for further projections: 

(5.9) = $F=l. 

Proposition 5.1. Under the given assumptions: (i) for solutions in Theorem ^.i[ i), (15.61) 

holds with oo(t) > for r>l, and then the matching with the boundary layer behaviour 
in (14.21) requires 

(5.10) iw^ 1 as r^+oo =S- p(s)=a (r)(l + o(l)). 

(ii) In particular, these are true for nontrivial nonnegative solutions. 

Proof, (i) follows from the construction of the boundary layer, 
(ii) This follows from Theorem I4.1( ii). □ 

Thus, projecting the PDE (15. 4ft onto the centre subspace of B* (i.e., by multiplying in 
L 2 by ipo(y) = F(y)) yields, for the leading mode ao(r), the following "ODE": 

(5-11) a'o = v y \ y=(p{T) ipo{ ( P( T )) + (K(v)v,if> ). 

The convergence (15.71) . which by a standard parabolic regularity is also true for the spatial 
derivatives, yields, as r — > +oo, 

(5.12) v y \ y=tp{T) = p( S )^(r) 7 i(l + o(l)) = a (rMr) 7l (l + o(l)), 7l = ^(0) = \. 

Finally, we need to estimate the last term in (15. lip : by (15. 6p . using that /t(a (r)) ^ 
for any ao(r) ^ via (ll.5p . there holds 

(5.13) (k{v)v, ipo) = («(a )a , F)(l + o(l)) = n(a )a (l + o(l)). 
Indeed, since J F = 1 for the Gaussian (1A. 13[) . in the last estimate, we have 

ip OO 

fF(y)dy=l- jF(y)dy=l-Oae-* 2 / 4 ) as <p = y{r) -> +oo. 

ip 

Thus, bearing in mind all above assumptions and estimates for generic patterns includ- 
ing (I5.6p . (15.101) . (15. 7p . and (I5.13p . we obtain the following asymptotic ODE for the first 
expansion coefficient cto(r) ^ 0: as r — > +oo, 



(5.14) 



ao(r) 40r 

IT 



h- y{r) + o(l)) + «(ao(r))(l + o(l)) . 



One can see that, by assumptions (11.51) . all the solutions of the non-autonomous ODE 
(15.141) are well defined for r 6 [0,+oo). Moreover, by classic comparison/monotonicity 
results for ODEs (S.A. Chaplygin's theorem of 1920s [H)]), it follows that, under the above 
hypotheses, solutions of ( 15.141) satisfy: 

(5.15) a (0) ^ ao(r) ^ for all r > 0. 
Therefore, we can always consider positive orbits: 

(5.16) oo(t) > for all r > 0. 

This makes our further asymptotic analysis easier. In particular, in view of ( 15. 15ft and 
( II. 5p . we can always omit all higher-order terms appeared via the above asymptotics. 

5.2. ODE regularity criterion. It follows from (15.141) . that a natural way to formulate 
a regularity criterion for the parabolic PDE (13. ip is to use the "ODE language"!! 

Theorem 5.1. (ODE regularity criterion) In the parabolic problem (13. ip . the origin 
(0,0) is regular, iff is globally asymptotically stable for the ODE (15.141) . i.e., any 
solution of (15.141) is global and satisfies 

(5.17) a (r) — )• as r — > +oo, i.e., In |a (r)| — > — oo. 

Proof, (i) Necessity. Given any classic solution u(x,t) (13.11) . one can always construct 
positive and negative Earner solutions u±(x,t) such that 

(5.18) U-(x, t) < u(x, t) < u + (x, t) in Q 

by standard comparison (Maximum Principle) arguments, [16]. Since, by Theorem I4.1( ii) 
and Proposition 15.11 such non sign-changing solutions u±(x,t) do obey our matched 
asymptotics, their positive (resp., negative) first Fourier coefficients satisfy the asymptotic 
ODE ( 15.141) for r 3> 1. Hence, by the BL-construction, (I5.17P implies that u±(x,t) — > 
as t — > 0~ uniformly, so, by comparison (I5.18p . the same does an arbitrary u(x,t). 

(ii) Sufficiency by contradiction. Let there exist a solution {a (r)} of (I5.17P (by (I5.15p . 
we may assume it to be positive) such that 

(5.19) lim sup^.,.^ ao(r) > 0. 

Then, by the ODE comparison, the same is true for solutions of (I5.17P with arbitrarily 
large Cauchy data at r = 0, i.e., for any 

(5.20) a (0) > a (0). 

Therefore, there exists a sufficiently large positive solution u + (x,t) of (13. ip . whose first 
Fourier coefficient satisfies ( I5.14p and (I5.20p . so the regularity is violated by (15.191) . □ 



In fact, this is quite natural and unavoidable: for semilinear PDEs, characteristic point regularity 
depends on asymptotic properties of ODEs, i.e., regularity issues for infinite-dimensional dynamical 
systems are characterized by ID ones. This reveals a sufficient and successful reduction of dimensions. 
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For the heat equation (ll.3p . with n = 0, integrating (I5.14p immediately yields 



(5.21; 



00 2, \ 



K 



: (0, 0) is regular iff f if(r) e 4 dr = +00 



which is indeed another equivalent form of Petrovskii's criterion (12. 7j) (in Khinchin's form). 

5.3. Applications: further regularity results. We now present a few corollaries of 
Theorem 15. 11 with simpler and more traditional conditions of regularity/irregularity. 

First of all, it follows from the ODE (I5.14p (and actually is true by comparison) that 
negative coefficients k(v) can "improve" the regularity of (0,0). Moreover, in this simpler 
case, we find a condition, under which any backward parabola has a regular vertex. 

Proposition 5.2. Let k{u) satisfy (ll.5p and let 

(5.22) k(u) < for u > 0. 

Then, for any backward parabola dQo with arbitrary tp's in (13.21) . (13.31) . the vertex (0,0) 
is regular. 

Proof. It follows from (15.141) that, for r> 1, 

(5.23) ^<_| K ( ao )|(i + (i))<_i| K ( ao )|. 

Then, on integration, assuming, without loss of generality, that ao(0) = 1, and checking 
an Osgood-Dini-type condition 

(5-24) / TOT = °°' 

0+ 

which obviously holds for the coefficients (II. 5p . we have 

(5.25) / ^|l_>i^ +00 as T ^oo. 

ao(-r) 

Hence, (15.241) reinforces (15.1 Tf) to hold. □ 

Secondly, for positive coefficients k, regularity can be destroyed. We first state the 
result establishing the conditions on monotone n(v ) > 0, under which the nonlinear term 
changes regularity for the pure heat equation into the irregularity. 

Proposition 5.3. Let n(u) satisfy (11.51) and let 

(5.26) k(u) > be increasing for u > 0. 

Let (15.211) be valid, i.e., (0,0) is regular for the heat equation (II .3p for N = 1. Denote by 
clq^t) — > as t —7- 00 the corresponding Fourier coefficient satisfying (15.141) for k = 0; 

(5.27) a (T) = a (0)e-^ / > (s)e ~" 2(3)/4ds for r » 1. 
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Then the linear regularity criterion (15.211) fails for the semilinear problem (13. ip and (0, 0) 
becomes irregular provided that the nonlinearity k is such that 

+00 2, > 

(5.28) / [ " i^F <P(r)e—f- + «(a (r))] dr > -oo. 

Proof. One can see that, in the present proof of a sharp estimate, one can omit both 
o(l)-terms in (15. 14f) . meaning that one can replace those by 1 + e and 1 — £ with an e > 
respectively and pass to the limit e — > + . 

As the first iteration of the full ODE (I5.14p . we have, for r 3> 1, 

(5-29) 4>__^( r ) e -^ =j> a (r)>a (r). 

In view of (I5.26p . (15.281) . we then obtain via the second iteration of (15.141) : 

(5.30) J>__J ;rV ,( T )e- 3 ^ + «(a (r)) for r » 1. 

Integrating this yields, by (I5.14p . that (0,0) is no more regular for such nonlinear coeffi- 
cients k(v). □ 

Corollary 5.1. Under the conditions of Proposition [37^ Petrovskii's backward parabola 
with the 2i/log log-factor H3.6I) is no more a regular vertex ofQo for the semilinear problem 
(13. ip provided that 

(5.31) k(v) > |ln|v||5 e - (3v ^ |ln|, ' ll)2/3 as f^0. 

Thus, (I5.3ip is the estimate, where the function in (II. 6p comes from. 
Proof. It follows from (I5.14p with n = that the function (15.271) reads 

(5.32) a (r) ~e"^ (lnT)3/2 as r ^ oo. 

Substituting (15.321) into (15.281) and changing the variable a (r) = v yields (15.3 1 p . □ 

Further iterating inequalities such as (I5.30p . one can obtain other sufficient conditions 
of the origin irregularity. For instance, if the integral in (I5.28P still diverges to — oo, 
integrating (I5.30p gives the next iteration estimate: for r 3> 1 

(5.33) a o (T)>a { i\r)^ao(0)eo l ^ K Ji \ 

where C\ > is some constant. Then, the next iteration leads to an awkward looking 
inequality: 

(5 ' 34) x / /[-^^^(ar**^™*)]*, 

+k a (0)e° 



Integrating it gives an estimate of a (r) > a (r) for r ^ 1 from below to be used also 
for the purpose of the irregularity, if a^\r) -ft as r — > +oo. If this fails, we then apply 
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the third iteration of the ODE (I5.14p again leading to a sharper estimate from below for 
the regularity, etc. 

Since the number of such iterations can increase without bound (and hence the same do 
the numbers of exponents and corresponding integrals in the argument of k(-) in (15. 34ft . 
etc.), it seems inevitable that a single and a simply finite integral criterion of irregularity, 
similar to the Petrovskii one (I5.2ip , cannot be derived for the nonlinear dynamical system 
(I5.14p in the maximal generality. In other words, the ODE criterion of Theorem 15.11 is a 
right way to regularity issues and is even optimal. 

In more general cases of equations in (13. ip . where, in our notations, 
(5.35) k = n(x,t,u,u x ), 

the derivation of matched asymptotics remains the same. The only difference is that, in 
accurate estimating of the integral in the last nonlinear term in (15. lip , we should take into 
account that v y ~ in the whole inner region due to the "centre subspace expansion" 
(15. 6p . so actually we integrate there k(-,0). But this term must also include integrals 
over the boundary layers close to y = ±y?(r), where the solution v and its derivative v y 
is sharply given by (15. 7p with the matching condition (I5.10p . We do not perform these 
general and, at the same time, rather straightforward and not that principal computations 
here, and restrict our attention to a particular model: 

5.4. Equations with a gradient-dependent nonlinearity. Let us very briefly con- 
sider equation (II .7p . Then, the first rescaling (13 .7p gives the equation 



\2 



V) 



(5.36) v T = B*u + k(v)v(v. 

It is easy to check that the BL-analysis yields the same asymptotics as in (15 .7p . with a 
similar proof. However, the eventual derivation of the ID dynamical system for the first 
Fourier coefficient ao(r) is now different: the nonlinear term is much weaker, since v y ~ 
on the centre subspace patterns, except a (^y) -neighbourhood of the boundary point 
y = (f(r). Overall, the nonlinear perturbation in (I5.14p is estimated as follows: 

J(oo) = (K(ao)agfe(^(l- f))] 2 (-^ 2 , Mv)) 



(5.37) 



K(a )al<p^f[ g ' (^(l-^)] 2 e-y^dy 
o 



where i/jq = F given by fl A. 13[) . Using the BL-profile (14. lip and setting z — - yields 



i 



J(a ) = K(a )afo 3 O- / e ->?fi-) e -W. dz 

o 

= K (a )a^ 3 ^e-^/e^f)dz. 



(5.38^ 



o 

Estimating roughly the last integral as follows: 

fe^'d-Vdz < eV, 
o 
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we obtain the following approximate dynamical system for ao( r ) > 0: 

(5.39) t<-^ <p(r) e'^ 4 + ^ ^0)^3(^-^/4 + ... . 

This is enough for us to prove that the nonlinear perturbation is now much weaker than 
that in (1514] : 

Proposition 5.4. For (I5.39p . Petrovskii's double log-function ( 13. 6 p forms a regular vertex 
(0,0) /or any function k(u) satisfying (11. 5p . 

Proof. Assuming that the linear term is dominant that creates the behaviour (I5.32p . one 
can check that, on this ciq(t), the nonlinear term in (I5.39P is always negligible, so (I5.30P 
follows. □ 

5.5. Backward paraboloid in R . More carefully, aspects of checking regularity of the 
vertex of a backward paraboloid in M. N was done in [20], where the authors applied match- 
ing techniques to the Navier-Stokes equations in M 3 . Now we present a few comments. 

For the iV-dimensional case ( II. 8p . the lateral boundary of the domain Qq in R N+1 is 
given by a backward paraboloid of the form 

(5.40) dQ : yjYlti = <f(r), r = -ln(-t), a, > 0, Ea? = l. 
Then a boundary layer close to the rescaled (via (14.11) ) boundary 

(5.41) dQ : E«ik| 2 = ! 5 

leads to a linear elliptic problem, which can be solved. Moreover, in the direction of the 
unit inward normal n to dQo, the boundary layer profile go(0 remains one-dimensional 
depending on the single variable 

(5.42) n = £ ■ n, 

so that g = go(rj) is still given by (14. lip . Therefore, in the expanding domain with the 
boundary 

(5.43) dQo(r) : Yl a i\Ui\ 2 = ( P( T ) ~ > +°° as t — )• +oo, 
the BL-profile is expressed in terms of the distance function: 

(5.44) 9q(v,t) = 1 - e -^W dist {^Qo(r)}. 

This allows us to apply the same blow-up scaling and matching techniques. 

The final ODE for ao(r) takes a similar to ( I5.14p form, with (p in the first term replaced 
by ip N , in a full accordance to ( I2.8p . However, the computations get more involved 
and further coefficients of this asymptotic ODE will essentially depend on the geometric 
shape of the backward paraboloid (I5.40p in a neighbourhood of its characteristic vertex 
(0, 0). However, final regularity conclusions remain approximately the same as for N = 1, 
including both cases of nonlinearities in ( II. 8p . 
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Appendix A. Hermitian spectral theory for operator pair {B, B*} 

For the maximal generality and further applications, we describe the necessary spectral prop- 
erties of the linear 2mth-order differential operator in M. N 

(A.l) B* = (-ir+ 1 A--^yV J/ , 

and of its adjoint B in the standard L 2 -metric given by 

(A.2) B = (-l) m+1 A™ + ^ y ■ V y + ^ / (/ denotes the identity). 
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Both operators occur after global and blow-up scaling respectively of solutions of the poly- 
harmonic equation 

(A.3) ut = -(-A) m u in 1^ x R+. 

Of course, for m = 1, (jA.lh and (IA.2j) are classic Hermite self-adjoint operators with completely 
known spectral properties; see e.g., Birman-Solomjak [9\ pp. 44-48]. However, for any m > 
2, both operators (jA.ip and (|A.2p . though looking very similar to those for m = 1, are no£ 
symmetric and do not admit a self-adjoint extension, so we follow more recent paper [11] in 
presenting necessary spectral results. In what follows, we mainly must concentrate on the less 
known case m > 2, naturally assuming that, for the classic self-adjoint case m = 1, we can 
borrow any result from several textbooks and/or monographs. 

A.l. Fundamental solution, rescaled kernel, and first estimates. We begin with the 
necessary fundamental solution b(x,t) of the linear poly-harmonic parabolic equation (|A,3p . 
which is of standard similarity form and satisfies, in the sense of bounded measures: 

jV 1 , 

(A.4) b(x,t) = t~2^F(y), y = x/t^ such that b(x, + ) = 5(x), 

where 5(x) is Dirac's delta. The rescaled kernel F = F(\y\) is then the unique radially symmetric 
solution of the elliptic equation with the operator (]A.2j) . i.e., 

(A.5) BF = -(-A) m F+^yVF + ^F = inR N , with / F = 1. 

For m = 1, F is the classic positive Gaussian 
(A.6) F (y) = _i F7ie -M 2 /4 >0 in R tf 

For any m > 2, the rescaled kernel function F(|y|) is oscillatory as \y\ — > oo and satisfies the 
estimate (for m = 1, this is trivial, with a = 2) [12\ I15j 

(A.7) \F(y)\ < D e~ d °\y\ a in R N , where a = G (1,2), 

for some positive constants D and do depending on m and N. 

A. 2. Sharp estimates in one dimension. For further use in our regularity study, we need 
some sharp estimates of the rescaled kernel, which we present for N = 1, where the regularity 
analysis gets also rather involved. Taking the Fourier transform in (|A.5P leads to the expression 

oo 

(A.8) F(y) =a J e~ s m cos(sy) ds, 

o 

where «o > is the normalization constant, and, more precisely |13j . 

oo 

(A.9) F{y) = 1 / e" s2m ^\y~\ J_i(s\y\) ds in M, 

o 2 

where J;, denotes Bessel's function. The rescaled kernel F{y) satisfies (|A.7p . where do admits 
an explicit expression; see below. Such optimal exponential estimates of the fundamental solu- 
tions of higher-order parabolic equations are well-known and were first obtained by Evgrafov- 
Postnikov (1970) and Tintarev (1982); see Barbatis [71 [8] for key references and results. 
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As a crucial issue for the further boundary point regularity study, we will need a sharper, 
than given by (1A.7|) . asymptotic behaviour of the rescaled kernel F{y) as y — > +00. To get that, 
we re-write the equation (|A.5P on integration once as 

(A.10) (_l)n»+ljp(2m-l) + _L y F = in R. 

Using standard classic WKBJ asymptotics, we substitute into (|A.10|) the function 

(A.ll) F{y) ~ y- 8 ° e ay " as y -»• +00, 

exhibiting two scales. This gives the algebraic equation for a, 

(A.12) (-1T M 2 — 1 = ^ and «5 = > . 

Note that the slow algebraically decaying factor y~ s ° in (jA.lip is available for any m > 2. 
For m = 1, this algebraic factor is absent for the exponential positive Gaussian profile 

(A.13) F ( 1 /) = ^-e-« a / 4 (m = iV = l). 

By construction, one needs to get the root a of (|A.12p with the maximal Re a < 0. This yields 
(see e.g., [8] and [HI p. 141]) 

( A - 14 ) «=^[-™(3^)+i«»(ap^iy)]=-* + i6o (db>0). 

Finally, this gives the following double-scale asymptotic of the kernel: 

(A.15) F(y) = y- So e~ d ° ya [C x sm(b y a ) + C 2 cos(b y a )] + ... as y -»• +00, 

where Ci^ are real constants, |Ci| + | C2 1 7^ 0. In (1A.15|) . we present the first two leading terms 
from the m-dimensional bundle of exponentially decaying asymptotics. 

In particular, for the linear bi-harmonic operator in (IB.lj) (N = 1), we have 

(A. 16) m = 2: a=|, do = 3-2"T, 6o = 3§-2 - T, and £0 = 3- 

A. 3. The discrete real spectrum and eigenfunctions of B. Both linear operators B and 
the corresponding adjoint operator B* should be considered in weighted L 2 -spaces with the 
weight functions induced by the exponential estimate of the rescaled kernel ()A.7p . We again 
more concentrate on the non-self-adjoint case m > 2, and refer to [9] for the classic one m = 1. 

For m > 2, we consider B in the weighted space L 2 (R ) with the exponentially growing 
weight function 



(A.17) p(y) = e a ^ a > in 



where a G (0, 2do) is a fixed constant. We next introduce a standard Hilbert (a weighted Sobolev) 
space of functions H$ m (R N ) with the inner product and the induced norm 

(v,w)p = I P(v) £ D k y v{y)D^w{y) dy, \\v\\ 2 p = f p(y) £ |^(2/)| 2 dy. 

IRJV fc=0 RiV k=0 

Then iJ 2m (R JV ) C L^(R^) C L 2 (R iV ), and B is a bounded linear operator from ff 2w (R JV ) to 

L^R^). The necessary spectral properties of the operator B are as follows [11] : 
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Lemma A.l. (i) The spectrum o/B comprises real simple eigenvalues only, 
(A.18) <r(B) = {A/3 = ~2^, fc = |0| =0,1,2,...}. 

(ii) The eigenf unctions ifip(y) are given by 

(A.19) Mv) = DPF(y), for any |/3| = k 

and form a complete set in L 2 (R) and in L 2 (R). 

(iii) The resolvent (B — A/) -1 for A c(B) is a compact integral operator in L pK 

By Lemma [A. 11 the centre and stable subspaces of B are given by E c = Span{V>o = F}, 
E s = Span{?/>/3, |/3 1 > 0}. Note also that the operators B has the zero Morse index, i.e., no 
eigenvalues have positive real part. In the classic Hermite case m = 1 (the only self-adjoint 
case), the spectrum is again given by (|A.18j) and the eigenfunction formula (IA.19P with the 
rescaled kernel (|A. 13|) generates standard Hermite polynomials; see [21 p. 48] for a full spectral 
account for the operator B. 

A. 4. The polynomial eigenfunctions of the operator B*. We now consider the adjoint 
operator (lA.ip in the weighted space L 2 *(R N ) ((•, -) p * and || • || p * are the inner product and the 
norm) with the "adjoint" exponentially decaying weight function 

(A.20) p*(y) = i =e-\y\ a > 0. 



We ascribe to B* the domain Hp™(R N ), which is dense in L^»(R Ar ), and then 

B* : H^(R N ) L 2 p *(R N ) 

is a bounded linear operator. B is adjoint to B* in the usual sense: denoting by (•, •) the inner 
product in the dual space L 2 (R N ), we have 

(A.21) (Bv,w) = (v,B*w) for any v € Hp m (R N ) and w G H^(R N ). 

The eigenfunctions of B* take a particularly simple polynomial form and are as follows: 

Lemma A. 2. (i) cr(B*) = cr(B). 

(ii) The eigenfunctions ip*p(y) o/B* are polynomials in y of the degree |/3| given by 

(A.22) r,(y) = ^t[y p + E [ jT m] jr( Ai»<V 

and form a complete subset in L 2 *(R N ). 

(iii) B* has a compact resolvent (B* - A/)" 1 in L 2 p *(R N ) for A cr(B*). 

Of course, for m = 1, (|A.22p yields standard Hermite polynomials, so, for m > 2, we call 
(|A.22p generalized Hermite polynomials. The bi-orthonormality condition holds: 

(A.23) (v^;) = <v 

Remark on closure. This is an important issue for using eigenfunction expansions of solutions. 
Firstly, as is well-known, for m = 1, the sets of eigenfunctions are complete and closed in the 
corresponding spaces, [9]. 
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Secondly, for m > 2, one needs some extra speculations. Namely, using (|A.23p . we can 
introduce the subspaces of eigenfunction expansions and begin with the operator B. We denote 
by L 2 the subspace of eigenfunction expansions v = ^cpipp with coefficients cp = (v,ifip) 
defined as the closure of the finite sums {Y1\/3\<m cp'tpp} i n the norm of L 2 . Similarly, for the 
adjoint operator B*, we define the subspace L 2 * C L 2 *. Note that since the operators are not 
self-adjoint and the eigenfunction subsets are not orthonormal, in general, these subspaces can 
be different from L 2 and L 2 ,, and the equality is guaranteed in the self-adjoint case m = 1, 
a = \ only. For m > 2, in the above subspaces obtained via suitable closure, we can apply 
standard eigenfunction expansion techniques as in the classic self-adjoint case m = 1. 

For m = 2 and N = 1 (this simpler case will be treated in greater detail), the first "adjoint" 
generalized Hermite polynomial eigenfunctions are: 

My) = i, My) = y, My) = -75 y 2 , My) = \ y 3 , 

(A.24) , V2 7 

My) = ^4 (y 4 + 24 ) ; Mv) = ^m(v 5 + 12 °v)^ My) = ^ (y 6 + 3eo y 2 ), 

etc., with the corresponding eigenvalues 0, — |, — |, — |, —1, — |, — |, etc. 

Appendix B. Semilinear bi-harmonic equations 

B.l. Regularity problem setting. Here, we show how our approach can be extended to 
higher-order PDEs, e.g., for the semilinear bi-harmonic equations having similar nonlinearities, 
with also zero Dirichlet conditions on dQo and bounded initial data u$ in Qq (~){t = —1}: 

{k(u)u 
in Qo, u = Ux = on dQ , u(x, -1) = u (x). 
K(u)u(u x y 

Then, after a proper similar matching with a boundary layer, we again arrive a nonlinear dynam- 
ical system viewed as a "centre subspace" approximation of solutions in the space of generalized 
Hermite polynomials as eigenfunctions of a rescaled non-self-adjoint operator. We also discuss 
the regularity problems for backward paraboloids dQo in R N x [—1,0), where the IB VP reads 

( k{u)u 

(B.2) u t = -A 2 u + \ W in Qq, u = §| = on 9Q , 

I k(u)u\ V«| 

where n is the unit inward normal to the smooth boundary of the domain Qq n {t}. Further 
extensions to 2mth-order parabolic PDEs are also discussed. 

Thus, we now show that a similar sequence of mathematical transformations can be performed 
for the fourth-order semilinear bi-harmonic equations (jB.ip . 

B.2. IBVP. We again fix N = 1, i.e., consider (|B.1|) with the simplest nonlinearity (|1.4p . 
leading to the IBVP 

u t = -u xxxx + n(u)u in Q = {\x\ < R(t), -1 < t < 0}, 
(B.3) \ u = u x = at x = ±R(t), -1 < t < 0, 

u(x,0) = u (x) on [-R(-1),R(-1)], 
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where Uq(x) is bounded and satisfies uq = u' = at x = ±i?(— 1). 

B.3. Slow growing factor <p(t). Similar to (12.5p . we assume that 

(B.4) R(t) = (-t)i <p(r), where r = - ln(-i) -> +00 as £->-0~. 

Here, the main scaling factor (— i) 1 / 4 naturally comes from the bi-harmonic kernel variables (see 



(|3.7p and (|A.4|) ) . and 99 (r) > is again a slow growing function satisfying f|3.3j) . For "shrinking 
backward parabolae" with 

¥>(t), 9?'(t) — >• as r — > 00, 

the regularity in the linear case k = was proved by Mihailov in 1963 [49|. 150] : in a certain 
sense, this extended the Gevrey-like result (12. 2p for m = 1; see (12. 3D . 



B.4. First kernel scaling. By (|3.2p . we perform the similarity scaling 
(B.5) u(x,t) =v(y,r), where y = ( _^ l/4 . 

The rescaled function v(y,r) solves the rescaled IBVP 

V T =WV + k(v)v = -Vyyyy - \yVy + K(V)V Ul Qq = { | ?/ 1 < ^ (t ) , T > 0} , 

v = Vy = at y = ±<£>(t), t > 0, 
[ u(0, y) = v (y) = u (y) on [-fl(-l), J2(-l)]. 



(B.6) 



B.5. Boundary layer. Sufhciently close to the lateral boundary of Qq, we naturally introduce 
the variables 

(B.7) z = ^py, v(y,r) = w(z,t) => w T = w zzzz - \ zw z + ^ zw z + k(w )w. 
The BL-variables now read 

(B.8) £ = (p5(T)(l-z), (p%(r)dT = ds, and w(z, r) = p(s)g(£, s), 

where p(s) is a slow varying function, for which eventually (|5.10p will hold by matching. 
Substituting into ()B.7p yields the perturbed equation 



, N 9s = Ay - \ -^3 £y e - & (1 - -|^)^ 

(B.9) / 

-3^73^5- ^5 + ^/3 k(P9) 9, wh ere Ay = -gM + \ g' . 

In this boundary layer, we are looking for a generic pattern of the behaviour described by (|B.9P 
on compact subsets near the lateral boundary, 

(B.10) |£| = o(y~3(r)) =>• |z - 1| = o(^"l (r)) as r ->■ +00. 

We next pose the same asymptotic behaviour (|4.8|) at infinity. Assuming that, by (|B.8[) . the 
rescaled orbit {y(s), s > 0} is uniformly bounded, by parabolic theory [12], we can again pass to 
the limit in (]B.9P in the asymptotically small perturbations, along a subsequence {sk} — > +00. 
Therefore, uniformly on compact subsets defined in (|B.10p . as k — )■ 00, 

(B.ll) y(s& + s) — )• /i(s), where ft s = A/i, /i = hg = at £ = 0, /i|^ =+OC) = 1. 
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Boundary Layer, m=2: the first stationary pattern g {%) 



Bo© 1 




FIGURE 1. [18] The unique stationary solution go (£) of the problem (|B.13p . 
The limit equation obtained from (|B,9|) . 



(B.12) 



h. 



Ah = -h^tt + \h, 



is again a standard linear parabolic PDE in the unbounded domain IR + , though now it is governed 
by a non self-adjoint operator A. Actually, we need to show that, in an appropriate weighted 
L 2 -space if necessary and under the hypothesis (|4.8p . the stabilization holds, i.e., the w-limit set 
of the orbit {h(s)} s> o consists of a single equilibrium: as s — > +oo, 



(B.13) 



(h(£, s) -»■ g (f), where Ag = for f > 0, 
I 50 = g'o = at £ = 0, 3o(+cx>) = 1. 
The characteristic equation for the linear operator A yields 
(B.14) _A 4 + IA = => Ai = and A 2 , 3 = ^ ( - \ ± i 

This gives the unique solution of (|B.13|) . shown in Figure [TJ 



(B.15) 



5o(0 = 1 



e 2 j 



[cos(l||) + ^sm(|g|)]. 



It turns out that the limit problem ()B. 12[) possesses a number of strong gradient and contrac- 
tivity properties. Namely setting by linearization 

(B.16) h(s) = go + w(s) =>• w s = Aw = —w^^ + \ W£, w = = at £ = 0, 
we arrive at the following (cf. Proposition 14.11 for m = 1): 

Proposition B.5. (i) (jB. 16|) is a gradient system in L 2 , and 

(ii) in the given class of solutions, the co-limit set f^o of (|B.16P consists of the origin only and 
it is uniformly stable. 
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Proof, (i) One can see that (|B. 16[) admits a monotone Lyapunov function obtained by multiply- 
ing W££ in L 2 : 

CB.17) U/W 2 = -/M 2 <0' 

Hence, (ii) also follows. □ 

Thus, quite similar to the second-order case, under given assumptions, we can pass to the limit 
s — > +00 along any sequence in the perturbed gradient system (|B.9p . Then, again similarly to 
m = 1, the uniform stability of the stationary point go in the limit autonomous system ()B.12j) in 
a suitable metric guarantees that the asymptotically small perturbations do not affect the omega- 
limit set; see [22^ Ch. 1]. However, at this moment, we cannot avoid the following convention, 
which for m = 2 is much more key than for m = 1, where the Maximum Principle makes this 
part of the analysis simpler, at least, for nonnegative or non-positive solutions (but for others 
of changing sign, this remains necessary). Actually, the convergence (|B.11|) and (|B.13|) for the 
perturbed dynamical system (1B.9|) should be considered as the main hypothesis, characterizing 
the class of generic patterns under consideration (and then (|4.8|) is its partial consequence). Since 
the positivity (negativity) is not an invariant property for bi-harmonic equations, a more clear 
characterization of this class of generic patterns is difficult. It seems that a correct language 
of doing this (in fact, for both cases m = 1 and m > 2) is to reinforce a "centre subspace 
behaviour" as in (|5.6p . rather than other (possibly, "stable") ones. Or, equivalently (and even 
more solidly mathematically), to impose the BL-behaviour (|B,13I) . 

Finally, we summarize these conclusions as follows: 

Proposition B.6. Under the given hypothesis and conditions, the problem (|B.9j) admits a family 
of solutions {called generic) satisfying (|B.13p . 

Such a definition of generic patterns looks rather non-constructive, which is unavoidable for 
higher-order parabolic PDEs without positivity and order-preserving features. However, we 
expect that (1B.13P occurs for "almost all" solutions. 

B.6. Inner region analysis: towards the dynamical system. As usual, in the Inner Re- 
gion, we treat the original rescaled problem (|B.6p . For simplicity of calculations, we again 
consider symmetric solutions defined for y > by assuming the symmetry at the origin: 

(B.18) Vy = Vyyy = at ?/ = 0. 

We next extend v(y,r) by for y > <p(r) and use the change (I5.2p . Since v = v y = on the 
lateral boundary {y = (p(r)}, one can check that, in the sense of distributions, 

V T = V T H, Vy = VyH, Vyy = VyyH, 

(B.19) Vyyy = VyyyH - Vyy\ y = ^5(y - If), 

Vyyyy = VyyyyH - v y yy\ y= ^5(y -if)- V yy \ y=(p 5'{y - 99). 
Therefore, v satisfies the following equation: 

(B.20) v T = B*v-Vyyy\ y=v 5(y-v)-Vyy\ y=v 5'{y-<p) + K(v)v in R + xR + . 
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Since such an extended solution orbit (|5.2p is uniformly bounded in L^*(R), we use the eigen- 
function expansion via the generalized Hermite polynomials (jA.22j) : 

(B.21) v(y,r) = Z (k > )^(r)r k (y)- 

Substituting (|R2ll) into pOuT) and using the bi-orthonormality property (1A.23P yields a dy- 
namical system: for k = 0, 1, 2, ... , 

(B.22) a' k = \ k a k ~ v yyy \ y=v{T) {5{y - <p(r)),ij) k ) - v yy \ y=tp{T) (S' (y - <p),if> k ) + ^h), 

where \ k = — | by (|A.18h . Here, < for all > 1. More importantly, the corresponding 
eigenfunctions V'fc(y) are unbounded polynomials and are not monotone for k > 1 according to 
()A.24j) , Therefore, regardless proper asymptotics given by (|B.22h . these inner patterns cannot 
be matched with the BL-behaviour such as (|4.8p . and demand other matching theory. However, 
these are not generic, so we skip them. 

Thus, we concentrate on the "maximal" first Fourier generic pattern associated with 

(B.23) k = 0: A = and if>^(y) = l (Mv) = F(v)), 

which corresponds to a "centre subspace behaviour" (15.6j) for the equation (lB.22p . which can be 
treated as another characterization of our class of generic patterns. The equation for ao(r) is: 

(B.24) a' = -Vyyy\ y=v{T) ipQ{if{T)) + v yy \ y=tp{r) %(cp(T)) + (As(ao)oo, V>o> + •• • 

Next, we use the boundary behaviour (|B.8p . ()B.13P for t> 1, which for convenience we state 
again: in the rescaled sense, on the given compact subsets, 

(B.25) v(y,r) = p{s)g Q {^{r){l - ^)) + ... , 

where qq is as in (|B.15p . Then, by the matching of both Regions for such generic patterns, (|5.10D 
must remain valid. Therefore, by (|B,25p . which by a standard parabolic regularity is also true 
for the spatial derivatives, we have that, as r — > +oo, 

2 2 4 

y=ip{j) ->■ p{s)<pz (r)7i -»■ ao(r)v?3 (r) 7 i, where 71 = ^'(0) = 2~s, 

y=lfi{r) -)■ -p(s)v?(r)72 -»• -a (T)(^(r)72, where 72 = ^"(0) = -i. 

Finally, for such generic patterns, we arrive at the asymptotic ODE for the first Fourier 
coefficient: 

(B.27) 



v yy\ 



(B.26) 

v yyy 



ao 



G±{ip{r),K) = 72V3(r)^o(^( r )) + 1W S (t)^'o( 1 p( t )) + K (ao) + ••• for r > 1, 



where, as usual, we omit higher-order terms relative all those remaining. Note that, for this 
ODE, the properties (|5.15p and (|5.16p remain valid by comparison. 

B.7. ODE regularity criterion and further applications. In general, the criterion of reg- 
ularity (I5.17P remains the same. However, it now reads: 

Theorem B.2. (ODE regularity criterion) In the fourth-order parabolic problem (IB.3p . the 
origin (0, 0) is regular in the class of generic solutions, iff any solution of the ODE (|B.27p 
satisfies (|5.17p . 
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Recall that by generic solutions we mean those that obey the boundary layer behaviour (|B.25p 
and hence, by matching with the inner region asymptotics, lead to the asymptotic ODE (|B,27j) , 
In this class, the proof of Theorem IB. 2 1 is straightforward. 

However, we must admit that we do not have a constructive way of describing generic solutions. 
In fact, this is not that exciting and/or surprising, since, even in the second-order case, solutions 
of constant sign were attributed to generic ones only by using the Maximum Principle, which is 
not available for bi-harmonic operators. For both second- and fourth-order parabolic equations, 
conditions of attributing solutions of changing sign to generic patterns are not fully known. 
Linear bi-harmonic equation: k = . However, the integrals in ()B.27[) are, in general, oscil- 
latory, so that a proper regularity analysis becomes not straightforward even in the linear case 
K = 0; see |18j . Then the following holds: 

oo 

(B.28) k = : J Gi{(p{T), 0) dr diverges to —oo <?=^ ao(r) — > as r — > +oo. 



Using asymptotic expansions of the kernel (|A.15|) and the corresponding eigenfunctions, as well 
as sharp values of the parameters (|A.16p . yields a more practical condition: 

(B.29) & = vHr)C 3 cos (bop! (T) + C 4 )e- do ^ /3 ^ + ... for r » 1, 

with some constants C^^ depending in an obvious way on C\ 2 in (|A.15|) and other parameters 
from (|A.16j) . Integrating yields 

f 

(B.30) ln|a (r)| = f tpl (s)d cos (b ipl (s) + C 2 )e- d ^ i/3 ^ ds + ... for r » 1. 

The regularity condition (|B.28|) is then re- formulated according to (|B.30P . Namely, the "critical" 
backward parabola occurs for the function (see |18[ § 7]) 

(B.31) ^(r) =3-3 2^ (lnr)2 + ... for r > 1, 

though, to guarantee divergence to minus infinity in (|B,28p . a special "oscillatory cut-off" of the 
function (r) may be necessary. 

Semilinear equations . For k / 0, instead of the linear (IB.29p . we deal with a nonlinear ODE 

(B.32) j£ =7^§(r)C 3 cos (& ^(t) + C 4 ) e"^ 4/3 M + K(a ) + ... for r » 1, 

and the analysis becomes more difficult. However, some of the results from Section [5.31 can be 
extended. Firstly, Proposition like 15.21 can be restored provided an oscillatory cut-off of v?(t) 
is performed for the first integral in the right-hand side of (|B.32p to be non-positive (though 
this business could look too artificial) . Secondly, a statement similar to Proposition 15.31 remains 
valid with the "linear" function (|5.27p replaced by 

(B.33) a (r) = a (0) e^ 3 («)«* M M+c 4 ) e -^ 3 M d Sj 

and with the corresponding changes in the integrals in (|5.27|) . (|5.28|) . 

Let us briefly (and more formally) derive "critical" nonlinearities k. It is easy to see that a 
somehow optimal and close to the critical dependence (|B.3ip is then achieved for the nonlinearity 

(B.34) K ( v) = _ n l^< for v^0+. 
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Indeed, solving the corresponding ODE without the linear term yields 

(B.35) a' = k(oo) / So(t ) j^y\ = r, where a (r) 0+ as r ^ oo. 

It follows from (|B.35P that, for the nonlinear coefficient (|B,34|) . 

(B.36) 5o(r) = e"^ 7 , 

so that, as is easy to see, that the linear term in negligible on the asymptotics ()B.36p . i.e., 
(B.37) yjl(r) e - d °^ /3 ^ = o(|k(5 (t))|) as r -> oo, provided that 

(B.38) 99(r)>(lnr)3 for r»l (cf. (IB^Tjl ). 

We thus arrive at a conclusion, which is similar to that in Proposition 15.21 for such negative n's, 
the vertexes of arbitrarily "wide" backward parabolae dQo are regular. 

Nevertheless, there are some principal differences with the much simpler second-order case. 
For instance, if the integral in (|B.30P diverges and both linear and nonlinear terms on the right- 
hand side of (1B.32P are sufficiently "balanced", i.e., both equally involved in the asymptotics of 
ao(r), the actual checking regularity/irregularity of the origin becomes a principally non-solvable 
problem. It is curious that the most interesting "interactional case" (of linear and nonlinear 
terms in (jB.32j> ) also begins at functions such as (II. 6|) . where the explicit constant 3y/n must be 
replaced by a more complicated one composed from those in (|A.16|) and 7^2 in (|B.26j) uniquely 
given by the BL-profile (|B.15|) , 

On the other hand, if the nonlinear term is asymptotically negligible on the "linear solutions" 
of (|B.32p . then the regularity and/or irregularity conditions remain practically the same as for 
the pure bi-harmonic flow. These are rather trivial results, which we do not intend to state and 
avoid such artificial "rigorous" theorems. 

B.8. Gradient dependent nonlinear perturbation. For the second equation in (jB.ip . the 

rescaled equation in (IB.6P takes the form 

(B.39) v T = B*v + k(v)v (v y ) 4 , 

so that using the same BL-profile (fBT27]) and the variables (jB~8|) . we obtain a similar dynamical 
system as in (lB.32p . where the weaker nonlinear perturbation is estimates as in Section [5.4} cf. 
(I5.37p . where tpo(y) = F(y) is the oscillatory kernel (|A.15p . One can complete these compu- 
tations; however, as before, such gradient dependent nonlinear terms do not affect the linear 
regularity criterion. 

B.9. Backward paraboloid in R . Again, in greater detail, regularity analysis in R^ for 
Burnett equations (with the bi-Laplacian rather than the pure Laplacian in the Navier-Stokes 
equations) is performed in [20 \ App. A], so we present here a brief notice only. For the equations 
(|B.2p . the lateral boundary of the domain Qq in M^"^ 1 can be given by the corresponding 
backward paraboloid 

(B.40) (E^i«iki| 2m )^ = (-i)^ fir), r = -ln(-i), a* > 0, £ af m = 1. 
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Again, a boundary layer study close to the rescaled (via (|4.ip ) boundary 
(B.41) 8Q : Zai\zi\ 2m = l, 

leads to a linear elliptic problem, which in the orthogonal direction becomes "quasi" one- 
dimensional, so that <7o(£) given in (|B.15|) depends on the single variable (|5.42p . Eventually, 
in Inner Region, the BL-behaviour leads to the stabilization to 

(B.42) g (y, r) = g (y,s (r) dist {y, 5Q (r)}) , 

where <7o(£) is as in ()B.15p . This makes it possible to derive the asymptotic dynamical system 
for the first Fourier coefficient and hence an ODE regularity criterion for generic solutions. 

The resulting asymptotic ODE for ao( r ) is similar to (|B,27p . with the extra multiplier ip N ~ 1 in 
the first two terms on the right-hand side. Inevitably, the final ODE will depend on the geometry 
of the backward paraboloid (|B.40p in a neighbourhood of its characteristic vertex (0, 0), which, 
in the most sensitive critical cases, makes it even less suitable for a definite regularity conclusion. 



B.10. More on generalizations. Using the above approach, there is no much principle dif- 
ferences and difficulties to treat the asymptotics of characteristic points for 2mth-order poly- 
harmonic equations 

(B.43) u t = {-l) m+l A m u + f{x,t,u,Vu,D 2 u,...) in Q , 

with zero Dirichlet (or others homogeneous) boundary conditions on OQq. Though, of course, 
some involved technicalities occur indeed. Since the first rescaled variables are 

(B.44) u(x,t) = v(y,r), y = { _ t) x 1/2m , t = -\n(-t), 

most interesting nonlinear terms in (|B.43P are now: 

/(") = (Zf) k(u)u, K(u)u\Vu\ 2m , K,(u)u\D 2 u\ m , etc. 

Then scalings (|B.44p lead to the rescaled parabolic equations 

k(v)v, 



(B.45) v T = B*v + < 



k(v)v \Vv | 2m , etc., 
k K,(v)v\D 2 v\ m , 

where B* is the linear adjoined operator (jA.ip . The corresponding dynamical systems for the 
expansion coefficients are obtained, as above, by (i) constructing a BL, and (ii) projecting 
the resulting PDEs ()B.45p . with the BL-approximation, onto generalized Hermite polynomials 
(|A.22p . The dynamical system, in general, becomes extremely oscillatory, and both linear and 
nonlinear terms can essentially affect regularity of the vertex (0,0). 
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Finally, we again mention that here our main goal: to show how the evolution of the first 
Fourier coefficient of generic solutions of bi-harmonic PDEs leads to an ODE regularity crite- 
rion, has been We must admit however that, in some cases, this did not end up with construc- 
tive/deterministic regularity conclusions, which are not always possible and are even illusive in 
general for higher-order nonlinear parabolic PDEs. 
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